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Abstract 



Versions of Tracy- Widom (TW) and Adler-Shiota-van Moerbeke (ASvM) approaches 
^ . are apphed to derive various partial differential equations (PDE) satisfied by joint eigen- 

K^ I value spacing probabilities of two coupled Gaussian Hermitian matrices (coupled GUE) . 

All the lowest (third) order PDE satisfied by the probabilities for the largest eigenvalues 
t~^ . of two coupled GUE are found, and the results of both approaches are compared. 

f^ I The TW approach allows to derive all PDE at once, while in the ASvM one starting 

t^ ' with different bilinear identities leads to different subsets of the full set of equations. 

An interesting result is that the joint probability for the largest eigenvalues of coupled 
Gaussian matrices satisfies a number of different PDE, and the previously known Adler- 
van Moerbeke equation (AvM) ^ is only one of them. Some of the new equations look 
like "coupled Painleve IV" and have usual Painleve IV equation as one-matrix limit, 
i.e. when the spectral endpoint of one of the matrices goes to infinity. This is in contrast 
C^ I to the AvM equation, which becomes trivial in this limit. Moreover, the new PDE, 

which stem from the matrix kernel approach of |23], do not contain derivatives w.r.t. 
the strength of coupling, unlike the AvM equation. In other words, they contain fewer 
independent variables and in this sense are simpler. 



* e-mail: igorrunianov@math.ucdavis.edu 



1 Introduction 

The problem of studying the joint eigenvalue probabihty distribution for coupled Gaussian 
matrices arises e.g. in studying Dyson process. A Dyson process is any process on ensembles 
of matrices in which the entries undergo diffusion. Consider ensemble oi n x n Hermitian 
Gaussian matrices. Then n curves describe changing eigenvalues with time. The probability 
that for each k the fc-th eigenvalue lies outside of a given set Jk at time t^ is given by the 
Fredholm determinant of a certain matrix kernel, called the extended Hermite kernel |23j . 
which is a particular case of a general matrix kernel first obtained by Eynard and Mehta [7j. 
Scaling this Gaussian process at the edge leads to the Airy process, introduced by Prahofer 
and Spohn [16J as the limiting stationary process for a polynuclear growth model. It is 
conjectured to be the limiting process for a wide class of random growth models, in particular, 
those belonging to the KPZ universality class [llj. Recently steps were made by Tracy and 
Widom [23] toward establishing it for the Asymmetric Simple Exclusion Process (ASEP), 
see e.g. [12], when they proved the long-standing conjecture that ASEP belongs to the KPZ 
universality class in a scaling limit of long time and large number of particles and its single 
point limiting distribution is the celebrated F2 [20]. In earlier works [221 123] Tracy and 
Widom (the approach in these works will be called TW further on) found systems of ordinary 
differential equations (ODE) for Airy, Gaussian and some other Dyson processes. Also, a single 
PDE for Gaussian process was known for a while from the work of Adler and van Moerbeke [3] 
(AvM), as well as PDE for Airy and Sine processes found from it in the corresponding scaling 
limits [3]. The last authors used two-dimensional Toda lattice (2-Toda) bilinear identity and 
Virasoro constraints to derive their equations. 

Here we compare the above two approaches to the joint distribution of the spectrum for 
coupled matrices. Even though we use a modified version of AvM approach, analogous to 
what we used for one- matrix ensembles in [T71 HB], which is closer to TW approach and is 
more flexible, we still find that system of equations derived from matrix kernel contains more 
information and allows in fact to obtain all possible PDE which such probabilities satisfy. 
This is quite similar to what we recently found for the general one-matrix case with several 
spectral endpoints [T9] . 

As Dyson first observed long ago [6], for Hermitian matrices with matrix elements inde- 
pendently executing Brownian motion subject to a restoring harmonic force, the equilibrium 
measure is the GUE measure of random matrix theory. With initial conditions at time ti 
distributed according to the GUE measure, the probability that at times t^ {k = 2, . . . ,m, 
we assume that ti < ^2 < • • ■ < tm) the matrix M(tk) is infinitesimally close to a value Mk is 
proportional to 

P ~ exp(-TrM,2) f[ exp f _IliM^^^h_lMt:ll'\ dM,--- dM^, (1) 

k=2 ^ ^ ^k-l J 



where c^ = e**~*'''+i. Alternatively, ([T]) can be interpreted as the equilibrium measure for a 
chain of m coupled n x n Hermitian matrices M^. The distribution of eigenvalues for this 
measure is expressible as the Fredholm determinant of an m x m matrix kernel |7] related to 
the kernel associated with the equilibrium random matrix ensemble - GUE in our case. It 
is derived by diagonalizing each Mk and then using the celebrated Harish-Chandra/Itzykson- 
Zuber (HCIZ) formula ^^: 

(where Xj, t/i are the eigenvalues of the Hermitian matrices X, Y, respectively, and, say, 
^(^) = nj<7(2^i ~ ^j) is the Vandermonde determinant), to integrate out the unitary parts. 
The result is the eigenvalue measure with a density P(xi, ■ ■ ■ ,Xr 



5 •^m) 1 
"^ ^ m— 1 



P ~ J]e VT^^+^;^'-^'='^ Yl det (e^^'^'^^'^^^M A{x,)A{xJ. (3) 

fc=l k=l ^ ^ 

It is shown in [71 |23] that for a chain of coupled matrices with probability density of this 
type the correlation functions can be expressed as block determinants whose entries are ma- 
trix kernels evaluated at various points, generalizing Dyson's expression for the correlation 
functions of a single matrix. In the case considered here the matrix kernel 

Kix,y) = iK,,ix,y))'^^^^ (4) 

is the extended Hermite kernel with entries 

f EfclJ e'=(*>-*^Vfc(a;)v5fc(y) if i > j, 
Kij{x,y) = <j (5) 

where (pk{x) = Pk{x)e~^ ^^ are the harmonic oscillator eigenfunctions and pk are the normal- 
ized Hermite polynomials. The actual kernel one is interested in is then 

Kij{x,y) = Kij{x,y)xjj{y), (6) 

where xj is the characteristic function of the set Jj, since the probability that for each k no 
eigenvalue lies in J^ at time tk is equal to det(/ — K^^). 

We restrict ourselves to the case of only two coupled matrices, m = 2, which we study in 
detail. Generalizations to the chain of several coupled matrices are quite feasible [TJ [5l 123] , 
but we defer them to a future work. 

Plan of the paper is as follows. In section 2 we show the normalization of the two-matrix 
integral and its one-matrix limit. Section [3] contains the description of the results. It connects 



the rest of the paper. In section HI using a version of the ASvM method [H E], we derive PDE 
for the joint distribution of two coupled Gaussian matrices from Toda lattice hierarchy. We 
use two copies of the AKNS system (see, e.g. |T5]) together with the 2-Toda (Liouville) 
equation and find PDE different from the one of [3], and also recover the AvM equation [3J. 
In contrast with it, our PDE have nontrivilal one-matrix limit. In section [5] we start from 
the matrix kernel and the line of TW [23], but then deflect from it to derive our version of 
first-order system of PDE of TW type, which employs matrix analogs of inner products u 
and w ^21] , shown in pTl [T8] to be universally related with the ratios of the corresponding 
1- Toda r-functions, i.e. one-matrix integrals of consecutive sizes. In the next section we derive 
matrix recursion relations for coupled case analogous to those in [18] for one-matrix unitary 
ensembles. In section [7] we transform the system of section \5\, find many first integrals and 
solve it, obtaining coupled analogs of Painleve IV equation [21] . In section [S] we match a 
certain combination of our TW-type equations from section [7] with a PDE from section [Hand 
find new correspondences among the quantities involved, in particular the relations between 
derivatives with respect to the coupling constant c arising from the Toda lattice approach 
([3] or section m here) and certain commutators involving the non-diagonal matrix elements 
of the matrix resolvent kernel [23]. In section Owe show how the system of Painleve IV- like 
PDE, the main result of section [TJ can be reduced to several PDE for the joint probability 
of largest eigenvalues wihtout auxiliary variables and with derivatives only w.r.t. the spectral 
endpoints of both matrices (i.e. without derivatives w.r.t. c). The last section is devoted to 
the conclusions. 

2 Matrix integral 

We study coupled Gaussian Hermitian unitary invariant ensemble, i.e. consider two-matrix 
integral^ 

Integrating out the "eigenvector" components of the matrices in the standard way, using 
Harish-Chandra/Itzykson-Zuber formula, and then transforming the integral, using the anti- 
symmetry properties of the HCIZ determinant, see e.g. [31 [23], one arrives at the formula for 
the joint eigenvalue distribution for this ensemble - the probability that all eigenvalues of the 
first matrix lie in a set Ji C M, while all eigenvalues of the second matrix lie in a set J2 C M: 



■ H N^ n(n— 1)/2 n 



n/ ^^.n/ %e-^"-?-^"^?+2^S"-^^'A(x)A(y), 

.•_i J J^ . — I J Jo. 



^' i=\JJ^ i=\J-J^ 



^Note that here, in the abstract two-matrix model, the normaUzation is different from the one for the 
Dyson process above. We will switch to the "physical" normalization later. 



where A(x) is the Vandermonde determinant over Xi. For brevity, we will denote r^^''^^ by 
just r;f. 

The coupling constant c varies in the range < c < 1, so that the integral converges and 
the physically identical domain — 1 < c < is excluded. 

First consider the case Ji = (—00,^1), J2 = (—00,00). Then change the variables i/i — )■ 
yi = yi- cXi to get 



/-, \ n{n-l)/2 n „g^ n ^^o 



A{y + ex) = YliKjiVi " Vj + c(a^j ~ ^j))- Due to the antisymmetry of the Vandermonde, the 
integral over new y variables decouples and turns into just a multiple Gaussian integral: 

Y\ dyie'^"y^A{y + cx) = c^^'^A{x)Y\_ dyie"^"^*. 



Therefore 



.•_i J —00 



=1 

Let 7 = 1 — c^, ij = 7X4, then 



-^K.) ~ n / ^^-^'' ■ rill = li^n / *.e-^;'.A=(x) 



.1/2 ryn(n^l)/2 ryn^/2 



i.e. such an integral is just a renormalized one-matrix largest eigenvalue probability: 

rii^i) = :^^n-"'^*"^(7ni). (9) 

For r-ratios of 2-matrix integrals of consecutive matrix sizes one gets therefore: 



Un{^^) = '-^ = -^.Ulr-'"'{^^l-i,l (10) 

W^n(ei) = %7^ = 7""'/'W^„'-'""*""(7'/'ei). (11) 



Thus, e.g. 



Fni^i) = U„Wni^i) = -F^™'^*"^(7^/2^i). (12) 

In the limit when ^i — )■ oo, we get the norniahzation for the defined above matrix integrals 
over the whole domain (see, e.g. [H] for the last formula in eq. (fT3|) below): 

^1— matrix 7r"/*^ n—i 
T- _ " l~matrix _ __^ TT 'i /-r on 

r2 -2(l-c2)- ^ ^ 

The above matrix integral is a r-function of 2-Toda integrable hierarchy [3] if one considers 
its modification by introducing two infinite sets of "times" or coupling parameters t and s, 

3 Description of the results 

We first apply a version of approach of [H [3] . We use 2-Toda bilinear identity, which supplies 
five simplest integrable PDE of its series - the 2-Toda (or Liouville) equation and two copies 
of AKNS system. Then the Virasoro constraints give us expressions of "time" derivatives in 
terms of spectral endpoint derivatives, which allow us to obtain our first system of five scalar 
PDE for joint gap probabilities of two coupled GUE w.r.t. the spectral endpoints: 

Theorem 1. The logarithmT of the joint spacing probability of two coupled Gaussian matrices 
or the joint two-time probability for the Gaussian Dyson process, together with the auxiliary 
functions U and W , satisfies the following system of PDE: 



AAT = 4c{UW -n/2), (E 

A^U = -2AoU - 2A^T ■ U, 1^ 

A^W = 2AoW - 2A^T ■ W, 1^ 

A^U = -2AoU - 2A^T ■ U, 
A^W = 2AoW - 2A^T ■ W. 
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Here we denoted 

A = B-i + cB-i, A = cB-i+B-u (15) 

Ao = Bo + c^Bo - V„ Ao = Bo + c^Bo - V„ V, = {1 - c^)cd. 



^CJ 



-1 ^ )5 Y^ „fc+l "^ 



iGendpoints of Mi iGendpoints of M2 



and 



' n n 

When all endpoints, say, of the second matrix go to infinity, the coupled matrix integral 
becomes the (renormalized) corresponding one-matrix integral, see the previous section. Our 
system of equations reduces in this limit to the three equations for single Gaussian matrices 
obtained in [IT]. The known AvM equation [3] for coupled Gaussian matrices, however, 
becomes trivial in this limit. So our system contains some additional information. 
Solving this system in general, we find two higher (fifth) order PDE satisfied by T: 

Theorem 2. The logarithm of the joint spacing probability for two coupled Gaussian matrices 
(or two time-point distribution for the Dyson process) satisfies the two higher-order PDE 
below: 

A (-A'F - AAoT - AA'TF + (■^^)^- ^ \ = _J (^^{A^TY + ^{Al + 2Ao)t] , (16) 



\2t? a a rr A A2rrT7 , y-^^) ^ \ _ A I ^ ^ A2rr\2 , 1/ J2 



A \ -A'F - 4AT - AA'TF + ^ ^- ]=~^\ ^(^^) + -(-^o + 2A)T 1 . (17) 

Expressions for F = UW , G = WAU - UAW, G = WAU - UAWG m terms of T read as 

I ~ n I ~ __1- 

F = —AAT + -, G = AT- —AAoT, G = AT - —AA^T, (18) 

Ac / Ic Ic 

so 



{AFf-G^ {A^ATY-A{AAqT -2cATY {AFf - G^ {AA^Tf ~ A{AA^T - 2cATY 



Ac{AAT + 2cn) F Ac{AAT + 2cn) 



Besides, we recover in our approach also the known third-order AvM equation [5]. 

Further on we consider in greater details the case of joint largest eigenvalue distribution for 
two coupled GUE. In this case, a simpler fourth-order PDE can be derived from the system of 
the first theorem, which we later compare with our equations derived from TW approach of 
matrix kernel Fredholm determinant. For this purpose, it is convenient to rewrite everything 
in a different way. Introduce notations: 

p+ = % + %, p_ = %-%, ^+ = ^1 + 6, e- = 6-6, (19) 

and 



(20) 



(21) 



1-c 

a = . 

1 + c 

The expression for F above now reads 

4(1 -a2) +2' 
Let 

G+ = WV+U - UV+W, G_ = WV^U - UV^W. 
The last quantities have expressions in terms of T, corresponding to flTSj) . 

G+ = 1/2V+T - l/4^+VlT - l/Ai^V+V^T + ^VJ)+T - l/2i+a^ ^^^~^-}^ , (22) 

Zc 1 — cr^ 



1 iV"^ -V'^)T 

G- = 1/2I?_T - 1/4:^^V^_T - l/A^+V+V^T Vj:)^T - l/2^_^^ -^. (23) 

2c (1 — o^) 

Then we have 

Theorem 3. The logarithm of the joint largest eigenvalue probability for two coupled Gaussian 
matrices (or two time-point distribution for the Dyson process) satisfies the fourth-order PDE: 

2FV+V_F - V+FV^F + G+G^ + 2F(^_G+ + ^+G_) + 8V+V^T ■F^ = 0, (EB 

where F is given by formula f21\) . while G+, G_ are given by formulas ^ME), ( fl^) - respectively, 
in terms o/ T = In r;f . 



Its one-matrix limit, i.e. the limit as ^i — )■ oo or ^2 ~^ 00, is a combination of the 3rd order 
equation in the derivative of T [211, which gives Painleve IV after integration, and Painleve 
IV itself. This is again in contrast to the equation f l82|) . trivial in this limit. 

It turns out that we can get more complete information from the matrix kernel approach 
of [23]. First, we derive the large system below, which turns out to be highly redundant. 
It contains 2x2 matrix variables analogous to the scalar variables of [2T]. Here the main 
variables are the entries of the matrix r, r^j = Rij{^i, C,j), which contains values of the matrix 
resolvent kernel R (for Fredholm matrix operator K"^) at the spectral endpoints, its diagonal 
entries are the first derivatives of the previously introduced T ~ In P, 

Tt = Trr = I^+T, r^ = Tr(o"3r) = V_T, 
while the anti-diagonal entries are expressed in terms of the matrix 



[^3, ^- l/2 [o-,r 
1-a 
where we have introduced the matrix a 



A = ^"^"rr.' ^ A = ia,a)A, (24) 



a = i-^as. (25) 

1 + c 

As a consequence of the system below, the scalar A^ {A^ ~ r^, square of the anti-diagonal 
part of r) is equal: 

2 _ V^n - V-rs _ Pilnr„^-PMnr„^ _ Adl^Jnr^j 

Theorem 4. The joint largest eigenvalue probabilities for the ensemble of two coupled Gaus- 
sian matrices and related auxiliary variables satisfy the following 2x2 matrix system of first 
order PDF: 

V+r = -^qP{I - a) - ^(/ + a)p| - K, r], mE) 

iasa)V^r = !#(/ - a) - ^{I + a)pl - [^, r], dUi 

V+q = -{^ + A)q + pU, ((Hi 

V_q = (-(73^ + A)q + aspU, m^ 



V+p = {^ + A)p-qW, (^T2^ 



V.p = (ag^ + A)p - asqW, m% 



^+1 = -q{^ + A) + up, dUi 

Vj=-I{a3^ + A)q + Upa3, ^B 

T^+P = k^ + A)- Wq, dUi 

V+U = Iq, (27) 

V_U = lasq, (28) 

V+W = -pp, (29) 

V^W = -pasp. (30) 

The other matrices in the above system of equations are degenerate (determinant zero). 

Besides, there are simple universal relations here which are direct analogs of the corre- 
sponding relations for the one-matrix case [18], involving matrix analogs u = Un = {<fn, {I — 
K^)ipn) and w = Wn = {'^n-i, {I — K:^)ipn-i) of inner products of [21], see sections [5] and [61 

Then, with G being the matrix with all elements equal unity (as in [23]) and Wn = 
(J - a)Wn{I + ff)/(l - a^), define 

f/ = _ 0n„e, ly = + QwnQ. 



Then in the Gaussian case U = y^n/2U and W = \/nj2W . 

Theorem 5. The correspondence between TW and ASvM variables found for 1-matrix case 
in IT1\ I7g| /. is modified for the coupled easel as follows: 



'i±lll^^ = det(/ - Qun) = 1 - Tr{eun), (31) 



^This is true in general for any matrix kernel of the form ((99)) like it was for the 1-matrix case [18 since 
we have not used the differentiation formulas for the harmonic oscillator wavefunctions here. 
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therefore 



^""'^^" ' = det(/ + Qijn) = 1 + rr(eiZ;„), (32) 



rr^ = 2%^, W = 2 '""Y/""\ (33) 

In fact, one has also a remarkable result: 

Lemma 1. 

Trifiw) = TrU ■ TrW. (34) 

All this plays a role in reducing the effective number of independent scalar variables and 
making the matrix system redundant. To solve it, we first reduced it to another matrix 
system (still highly redundant for the above reasons). Let 

X+ = qp + pq, X^ = qp- pq, (35) 

$ = 2{pUp - qWq), G = 2{pUp + qWq), (36) 

then 

Theorem 6. The following closed system of ten PDE is a direct consequence of the system 
of fourteen PDE obtained in section\^ 

V+r = -1/2X+ - l/4[or, X+] + l/4{or, X_} - K, r], ^M 

{<y^a)V^r = 1/2X_ - l/4{a, X+} + l/4[a, X_] - [e, r], ([123) 



and 



2)+X+ = <l>-[e + i,X„], (37) 

P+X_ = -[e + i,X+], (38) 

I?^$ = 3X2 -8nX+ + X2 +{^ + i,G'}, (39) 

I?+G' = {e + i,$} + [X+,X_], (40) 
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D_X+ = 1/2(^3,$} - [as^,X^] + [A,X+], (41) 

I?_X_ = l/2[a3, G] - [as^, X+] + [A, X_], (42) 

V_<^ = {as, Xl - 4r2X+ + X!} + X+aaX^ - X^agX^ + {^3^, G} + [A, $], (43) 

V^G = [0-3, {X+, X_} - 4r2X_] + X+(T3X_ - X^a^X^ + {^3^, $} + [A, G], (44) 

Proof. The first two equations above are just the first two equations of theorem HJ getting the 
others is straightforward expanding of their left-hand sides by definitions (13 5 p and (I36p and 
then applying the rest of equations of theorem H] together with definitions flM|) . 

D 

The system has a number of first integrals some of which are matrix while the other involve 
only some of the scalar components, thus forcing one to split the system into scalar parts 
eventually. As the main consequence of the system, we found a complete system of inde- 
pendent third-order PDE satisfied by joint largest eigenvalue distributions for two coupled 
GUE: 

Corollary 7. 



V, = V+XtV_Xt - 2XtX^F - GtG-s = 0, ([7^ 

D^ = Aa^iV+A^V^A^ - 2X3(^2)2) - A+A. = 0, ( [77^ 

2a^A^Vt = -2a'^A^V^ = -FDt = a^F^s, (45) 

FiV+A^A^ - V^A^A+) = A'^{G3V+Xt - GtV^Xt), ([17^ 



where 



Vt = {V+Xtf - F{2Xl + J) - Gl V3 = {V_Xtf - F{2Xl + J) - Gj, (46) 

Dt=AcT\V+Ay + 2a^A^J-Al, D3 = Aa^V^A^ + 2A^J - A^. (47) 

9.9 .9 (Pi -(721)2) In r^ 

Xt = -2V+rt - 2a^A^ = -2V^r3 - 2^2 = -2^^^ ^ -, 

1 — cr^ 
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J 

n 1 



X3 = -2D+r3 = V^rt = -2V+V^ Inr, 



D+ = V+A\ D_ = V_A\ 



Gt = Ht + A+, G, = Hs + A_, Ht = An-2i^V^n + i.X:,, H^ = Ar,-2^_V_r3 + ^+X3. 

A+ = tr{A, (X+) J, A_ = -tr{A, (X_) J, (48) 

({, } is an anti- commutator, subscript 'a' means the anti-diagonal part of a matrix). 

The last four expressions clearly have the three-term structure resembling that of Painleve 
IV equation for one- matrix Gaussian ensemble, as do the left-hand sides of fll73p and fll74p . 
And indeed the limit of Vt, V^ and Vx as either ^2 — ^ 00 or ^1 — )■ 00 is the Painleve IV 
equation itself. The one-matrix integrals for GUE largest eigenvalue probabilities, which are 
the appropriate solutions of the limiting Painleve IV equations, should in fact be considered 
as boundary conditions for the corresponding PDEs of the last corollary. 

This nice form of the equations is reached at the expense of still keeping two auxiliary 
variables A^ and A^ in them. They can be explicitly eliminated giving a smallest complete 
system of four independent third-order PDE for T, which, however, lack the nice Painleve 
IV-like structure of the previous ones, see section O 

Comparison of a combination of 4th-order PDEs of theorem [TD] and 3rd-order PDEs of 
corollary [7] obtained from the variant of TW system in section [7] with the 4th-order PDE of 
theorem [3] above from our version of AvM approach, gives the correspondences of the variables 
of both approaches, among them - some new ones, specific for the case of coupled matrices: 

cr3[r„,r'+rj = -2cdcrt, a3[r„,r'_rj = 2c(9cr3. (49) 

These last two new important relations show that the main additional variables appearing 
in the TW approach to coupled (Gaussian) matrices - the commutators on the left-hand 
side - directly correspond to the main new objects appearing in the ASvM approach to the 
same problem - the derivatives w.r.t. the coupling c (or w.r.t. time t, recall c = e~*) of the 
log-derivatives of the joint largest eigenvalue probability. 
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4 System of PDE from Toda lattice hierarchy 

Along with [3j, we first consider the following 2- Toda r- function: 



of section HJ From the bilinear identity for the 2- Toda hierarchy |25] , see also [3] , 



r„+i(t, s - [z])TUt', s' + [z])e^^^^^-^'''>-\'^-"^-'dz, 



2=0 



where [z] = {z, z'^/2, z^/3, ■ ■ ■) - infinite vector, for m = n ± 1 one obtains the series of PDE: 

7 = 2^ 7^ ^^^) 



dtkdsi ^ r, 

1=0 



and 



(9^ In r„ _ >;p Pii-ds)rn+i ■ Pk-i-i{ds)rn-i ,^^. 

1=0 

both of which give as the simplest {k = 1) equation the 2-Toda or Liouville equation in terms 



9tidsk "t^^ r, 

give as the e 
of r-functions: 



The identity for m = n gives another two series of PDE: 



(53) 



A; ~ ^ 

d Tn+l _ Y^ Pi{-dt)Tn ' Pk-i{dt)Tn+l , . 

dtk Tn j^ TnTn+1 



k 



--—In => , (55) 

which give the simplest nontrivial equations at k = 2, namely two copies of AKNS system: 
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dUr, d^U„ dHnT„^, dW„ d'^W„ S^lnr, 



ds2 dsf ds\ "' ds2 dsf dsf 

Then we use the Virasoro constraints [U El El [IT] , connecting the derivatives w.r.t. the spectral 
endpoints with the time derivatives of the integral r^it^s) ( l50l) . We need only the k = —1 
constraints: 

-2— +nh + J2 ItiQl 2c— \ T^{t, s), {51 

^ 1=2 ^ ^ ^ 

B.,Ti{t, s) = V_ir„^(t, s) = \2—-ns, + Y, l^i-. + 2c— \ rlit, s), (59) 



dsi -^-^ dsi-i dti 

^ 1=2 



and the k = constraints: 



BoT;i{t, s) = Voriit, s) = I -2 A + ^^^^ + E ^^'|; + ^^ | ^n{t, s), (60) 

BoT^it,s) = VoTiit,s) = I2A + ^^^^^ + f;/..^ + c|;| r;f(t,.). (61) 

So on locus tk = Sk = one gets: 



ds2 2 -^ dsi dc 



OSi Oti 

Constraints for the second derivatives follow from the above, using the commutativity of 
boundary and time derivatives. On the locus those we need reduce to: 
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Expressing the time derivatives in terms of the boundary derivatives as usual in the ASvM 
approach [H [3], we finally find for the second derivatives of Inr;^: 



a2_ln7^_ (g_i + cE.xfhiTi + 2(1 - c^)n 
dt\ 4(l-c2)2 ' 

a^lnr^^ _ (ci3_i + B^xfXTvTi + 2(1 - c^)n 
ds\ 4(l-c2)2 ' 

a^ lnr„^ (i3_i + ci3_i)(cS_i + ^_i) lnr„^ + 2c(l - c^)n 



(62) 
(63) 
(64) 



dt^dsx 4(1 - c2)2 

After plugging the last expression into the 2-Toda equation, the "boundary- Toda" equation 
for coupled matrices follows: 

(i3_i + cB^i){cB^i + B^i) lnr„^ = 4(1 - c^fUnWn - 2c(l - ^)n. (65) 

Taking the corresponding formulas for the ratios f/„, Wn-, and recalling the definitions from 
section [3], 

one obtains four "boundary" equations from two AKNS copies: 



A^Un = -2(1 - c2)2(i3o - cd,)Un - 2(1 - c2)((l + ^)n + ^)Un - 2A^ lnr„^f/„, (66) 

A^Wn = 2(1 - ^)\Bo - cd,)Wn - 2(1 - c2)((l + ^)n - l)Wn - 2A^ lnr„^iy„, (67) 

A^Un = -2(1 - c2)2(i3o - c9J?7„ - 2(1 - ^){(l + ^)n + c2)f/„ - 2^^ \^r^u,^^ (68) 

^2iy„ = 2(1 - ^)\Bo - c9e)iy„ - 2(1 - c2)((l + c2)n - l)Wn - 2A^ Inr'^W^. (69) 

Now let us pass to the "physical" Dyson process variables, namely substitute 

e ^ ^=, r„^ ^ c^^n , (70) 

V 1 — c^ 
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(if we did only the second substitution above we would get exactly the matrix integral of 
section 2 as r^) and so 

f/„ ^ c"f/„, Wn ^ Wnlc^-\ (71) 

which entails also 



d f d \ d \ ( d \ 



(72) 



anc?, importantly, change in the meaning of the partial derivative with respect to c (due to 
the acquired dependence on c of former ^ under the change of variables ( 1701) ). so that one has 
to substitute everywhere 



- cdc = 2c'd^ -^ 2c' d^ + r^i-^:^ . (73) 

After these transformations the previous system of PDE reads: 

^^lnr„^ = 2c(2(l-c2)f/„iy„-n), ([65]) 

A'Un = -2{Bo + c'BoWn + 2(1 - c2)c9,f/„ - 2AHnT;lUn - 2^{2n + 1)?7„, ([MD 

J^Wn = 2{Bo + c'Bo)Wn - 2(1 - c')cdcWn - 2A' In r^Wn - 2c\2n - l)Wn, (EZD 

A'Un = -2{Bo + c'Bo)Un + 2(1 - c')cd,Un - 2A' In r„^f/„ - 2c\2n + l)Un, (IS 

A'Wn = 2{Bo + c'Bo)Wn - 2(1 - c')cdcWn - 2AHnT^Wn - 2c2(2n - 1)W„. 

At last, it is convenient to remove the last terms in the last four equations by passing to new 
functions, 

U = Un{l - ^T^^l\ W = Wn{l - c2)-("-l/2), 

which also simplifies the coefficient in equation (1^3]) and the final system becomes that of 
theorem [H 

Now we try to solve this system. The procedure goes in parallel with one-matrix case 
(see [l8l[T9]). First introduce quantities 

F = UW, G = WAU - UAW, G = WAU - UAW, 
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Co = WAoU - UAoW, Go = WAqU - UAqW. 

Multiplying eqs. f l66|) and fl68|) by W and eqs. f l67|) and fl69|) by f/ and taking the appropriate 
combinations to express everything in terms of the newly introduced functions, we obtain our 
system in the form: 

AAT = 4c{F-n/2), (74) 

AG = -2AoF, (75) 

AG = -2AoF, (76) 



jAFy - G^ 
~2F 



j,^p _ y-^ ) ^ ^ _2G^ - AA'^T ■ F, (77) 



jAFf - & 
~2F 
We will need commutation relations among our operators: 



j,^F - ^ ^^ = -2Go - 4^'T • F. {n 



[A, Aq] = A, [A, Ao] = 2cA - c^A, (79) 

[i, Ao] = 2cA - c^A, [A, Ao] = A. (80) 

Equation (17^ is in fact just the definition of F in terms of T. Plugging it into the right-hand 
sides of (^^ and (175]) and applying the last relations allows one to integrate this first couple 
of linear PDE and obtain expressions for G and G in terms of T: 

AT ~ ~ AT 

G = -(A-1)^, G = -(A-1)^. (81) 

The last expressions already allow to recover the Adler-van Moerbeke equation [3l H] in our 
approach. It follows in fact from a trivial identity: 

AA\n- = AA\n-, 

rewritten as 

'G 




Plugging in expressions for F, G and G in terms of T, we get a single PDE for T = Inr;^: 



Theorem 8. The logarithm of the joint spacing probability for two coupled Gaussian matrices 
(or two time-point distribution for the Dyson process) satisfies the third-order PDF, 

/A^}}± = .l (-^°-l)^ (82) 

AAT + 2cn AAT + 2cn 

which is nothing but the equation first derived in [3] using two first {k = 1,2) bilinear identities 
of the Toda series, while we started from two couples of first identities of AKNS (non-linear 
Schrodinger) series instead of the second {k = 2) Toda identity. We, however, can derive also 
equations (fT6|) and (TT71) of theorem [2] from the system of theorem [H see Appendix A. It seems 
that in general these equations are hard to simplify. 

Further on we will consider the case of only one boundary point for each matrix, i.e. study 
the joint largest eigenvalue distribution. Let us then rewrite everything in a different way, 
which will be convenient for the comparison with the variant of TW approach below. 
Recalling the definitions of operators P+, P_, introduced in section [3l we get 

A=^-^iV^ + aV_), A=^-^iV^-aV_), 

A'~A^ = {1- c2)P+P_, A'' + A''= ^^^^{Vl + a^Vl), 

Bo + Bo = im+V+ + e_P_), B,-Bo = l/2(e+P- + ^-V^). 

The system we obtained in theorem [1] now acquires the following form, if we take differences 
and sums of its last four equations: 



{VI - a^Vl)T 



n 



m 



4(1 -cr2) 2' 

V+V_U = -{i+V_ + i-V+)U - 2V+V_T ■ U, (83) 

V+V^W = (^+I?_ + ^-V+)W - 2V+V^T ■ W, (84) 

{Vl + a'^Vl)U = -(1 + <y'^){i+V+ + i-V^)U + SacdJJ - 2(1?^ + a'^V'i)T ■ U, (85) 

{Vl + a^Vl)W = (1 + a^){^+V+ + ^^V^)W - SacdcW - 2{V\ + a^P^ )T ■ W. (86) 
Let 
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G+ = WV+U - UV+W, G- = WV^U - UV^W, G^ = WcdJJ - UcdcW. (87) 

Then from the first two equations for U, W we derive one linear and one non-hnear equation: 

I?+(G'_ + 2e_F)+I?_(G'+ + 2e+F)=0, (88) 

2FV+V^F - V+FV.F + (G_ + 2^_F)(G+ + 2^+F) = 4(^+^_ - 2V+V^T)F^. (89) 
From the longer pair of equations also one linear and one non-linear PDE result. They are 

I?+(G+ + 2e+F) + a'P_(G_ + 2e-F) = {l-a^){^+V+ + ^^V^)F + 8acdcF + 4{l + (T^)F, (90) 

2F{Vl + a^Vl_)F - {V+Ff - a\V_Ff+ 

+ (G+ + 2i+Ff + a\G^ + 2i^Ff - 2(1 - a^){i+{G+ + 2^+^) - i^{G^ + 2^^^))^ = 

= IQaFGc - 8{Vl + ^22)2) {-^t^F'^ + 4(^+^2 + C^)F'^. (91) 

The pair of linear equations here is, of course, satisfied by the previously found general 
expressions for G and G. Since now 

G=^i±^(G+ + aG_), G = ^^(G+-aG_), (92) 

they correspond to expressions [22] and 1231 for G^ and G-, given in section [3] before theorem 

El 

The last expressions, together with fl2T|) . are to be plugged into eqs. fl89|) and fl9T|) . After that 
eq. fl89|l becomes a single 4-th order PDE for T = Inr;^, that of theorem [31 

As for the equation f[9T|) . the obstacle to obtaining another PDE for T from it remains, 
since we have not found a convenient expression for Gc- We can get much more from the 
matrix kernel approach [23], which is the subject of the next sections. Nevertheless, the 
considerations presented here will be very helpful there to decipher the complicated system of 
matrix equations and to find more of the universal relations like those in [18] for one-matrix 
UE. 
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5 Matrix Kernel and TW equations for coupled Gaus- 
sian matrices 

We consider sets J^ = {C,k, C)o), i.e. the simplest case of only largest eigenvalues and defer the 
due generalizations to several spectral gaps for a future work. Then 

det(J - K') = P(M(ti) < 6, • • • , M(t„) < e^). 
We first follow [23] here and start with the general formula: 



d^K' = -KSk, (93) 

where Sk is the diagonal matrix with all entries zero except for the k-th and {5k)kk = ^{y — ^k}- 
Then we introduce the (matrix) resolvent operator R = K'^{I — K'^)~^ and see that 

9fclndet(J - K') = -Tr(/ - K'Y^dkK' = RU^u^ik)- (94) 

We keep most of the notation in [23| i.e. let D = d/dx, M = x- - the multiplication operator, 
p= {I - K-^)-^ = I + R,^ = diag(^fc), d^ = dmg{d^k), x{x) = diag(xj,(a;)), (5 = Efe 4, © - 
matrix with all elements equal unity; introduce matrices r, r^, Vy such that 

rij = Kij[t,i,t,jj, (,^xjij = [(yxK)ij[t,i,t,j), {^y)ij = yOyJ^jijyki^ki)- l9oj 

Define also (here our normalization is as e.g. in [21] for 1-matrix models rather than the one 



m ^ ^ 

^ = hH^^^n, ip = 6y_\(^„_i, (96) 

where 6„_i = \rn/2 for the Gaussian case, and introduce matrix functions 

Q = pip, P = pi), Q = (fXP, P = -ipXP, (97) 

and their values at the spectral endpoints 

Qij ^ij\<,i)i Qij ^ijK^jJi Pij ^ij\Si)-i Pij ijv^jl' 1'-'°/ 

Following [23], also consider instead of ([5]) a modified kernel with the same Fredholm deter- 
minant: 

f E::Je(^-"^^*--*^Vfc(a:)v^fc(y) if^>J, 
K,,{x,y) = l (99) 

I -Er=ne('""^^*'-*^Vfc(x)¥.fc(l/) if^<j. 
Using the well-known differentiation formulas for the harmonic oscillator eigenfunctions: 
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d'^k /—r difk-i i—r 

-xipk + V2fcv9fc_i, — - — = XLpk-i - v2kipk, 



dx dx 

we have 



{D + M)v9 = V2nV', {D - M)^ = -V2nip. 
Then the two operator identities called Lemma 3 on p. 24 of [23] read: 

{D + M)Kij - e''-'^Kij{D + M) = -2iP{x)^{y), 



-i^Z ^1 



{D - M)K,, - K,,{D -M) = -2^{x)i,{y). 



The last formulas multiplied on the right by the matrix x, together with relation [-D, x] =5 
(and [M, x\ = 0), if we let t = diag(tj) and so e* = diag(e*^), lead to 

e-'{D + M)K'^ - K^e-\D + M) = -2e''%l){x)Qx{y)^{y) + K5e-\ (100) 

e\D - M)K'^ - K-^e\D - M) = -2^{x)Qx{y)^{y)e' + K5e\ (101) 

Multiplying (IIOOI) and (110 II) by p on the left and on the right, replacing K'^ — )■ K^ — / on the 
left-hand sides and using the above definitions we get the two equations of Lemma 4 of [23] 
(with the correction of swapping the matrices e* and in the last formula below): 

e-\D + M)R - Re~\D + M) = -2P{x)e~'QQ{y) + R5e-'p, (102) 

e\D - M)R - Re\D - M) = -2Q{x)Qe'P{y) + R5e'p. (103) 

At this point we depart from the line of [23]. Recall now that we consider the case of 
two matrices, in other words only two time points. Therefore, denoting c = e"^*^"*^'* and 
introducing matrices 

(1 + c) ^ (1-c) (1 + c) , (1-c) 
eL = — ^ — /+ — ^ — CT3, eu = — ^ — / ^ — CT3, (104) 

where a^ is the usual 3rd Pauli matrix, we can rewrite (I102p and (I103P as, respectively, 

eL{D + M)R - ReiiD + M) = -2P{x)eLQQ{y) + RSclP, (105) 

eu{D - M)R - Reu{D - M) = -2Q{x)QeuP{y) + RSeup, (106) 

Take i,j entries in (llOSp and (I106P and set x = ^i, y = C,j. This gives 
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eiTx + TyCL = -eL^r + rei^ - 2peLQq + reir, (107) 

eu^x + ryCu = eu^r - reuC - '^qQeuP + reur, (108) 

Since ei + eu = (1 + c)/, ei — eu = {1 — cja^, adding and subtracting the equations (I107p . 
(llOSp . we get, respectively, 

Tx + ry-r'^ = -p{I + a)Qq - q&{I - a)p - [a^, r], (109) 

o'Tx + Tya — rar = —p{I + a)Qq + g9(/ — a)p — [^, r], (HO) 

where we have introduced the matrix a, 

a=^—^a;. (Ill) 

1+c ^ ^ ^ 

To derive our version of first order TW system, we introduce the matrix function m, 

u = {^x,Q) = {Q.^), (112) 

analogous to the scalar function u for one-matrix case. It has not been used by Tracy and 
Widom in [23] but had already appeared in their earlier work [22], where the matrix Airy 
kernel for the Airy process was considered. Similarly, we introduce the matrix w, the analog 
of the scalar function w from [2T] : 



w = {i,x,P) = {P,i')- (113) 

Acting by operator equation (llOSp upon function ip from the left, we get 



^AiD^M),p 



V = ^^"Y^P + ^)Q - vf ^(^ + ^) = -^(^ + ^)®^ + 1^^^^ + ^)^' 



I.e. 



(/ + cj){DQ + xQ) = 2P(/ + ^) ( y ^ - 0^) + R^i^ + ^)Q- (114) 

Similarly, acting by equation (I106p upon function ip from the left gives 

(J - a){DP - xP) = -2Q (J^ + 0(J - a)w{I - a)-^j (/ - (x) + R5{I - a)P. (115) 
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On the other hand, acting by fll06p upon ipx and by fllOSp upon ipx from the right gives, 
respectively. 



{DQ + Qy){I-a) = 2[ J^ - uQ ] {I - a)P + Q{I - a)5 ■ R, (116) 



and 



{DP - Py){I + a) = -2(/ + a) U^ + (/ + a)-'w{I + a)e j Q + P{I + cx)6 ■ R. (117) 



Next we use the fact following from equation ( 1931) : 

dkp = pdkK^p = -R6kP (118) 

to derive the formulas [23] : 

dP ^ dp ^ 

Similarly, we get the corresponding expression for ©.-derivatives: 

'D-Qij = ^3-i-i^i) - {rcr3q)ij, V^qij = —{^i)a3 - (gasr)^^, (121) 

dP ^ dp ^ 

V_pij = o"3-^(6) - {r(r3p)ij, V^pij = -j-i^i) - {P(^3r)ij. (122) 

We take x = C,i, y = ^j in the eqs. (I114p . (IllSp . (I116p and (I117p and use the previous formulas 
to obtain PDE with spectral endpoints ^1,^2 as indepenent variables: 

(/ + a)V+q = -(/ + a)^q + 2p{I + a) (^ - Qn) - [a, r]q, (123) 

(/ - a)V+p = {I- a)ip - 2q (^ + 6(1 - a)w{I - a)-'^ {I - a) + [a, r]p, (124) 

V+q{I -a) = -qi{I -a) + 2 Lj^ - uQ\ (/ - a)p - q[a, r], (125) 
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V^p{I + a)=pai + a)-2{I + a)U^ + {I + a)-'w{I + a)e]p + p[a,r]. (126) 

and the other four equations for V_ derivatives: 



(J + a)V^q = -(/ + a)a3^q + 2a3p{I + (t) U"^ - Bu ] + a-^ril + a)q - (/ + a)ra3q, 
i.e. 

(/ + a)V.q = -{! + a)a3^q + 2a;p{I + a) f W^ - Gm ) + [(T3, r]g, (127) 



and, similarly, 



(/ - a)V^p = (/ - a)a^iip - 2asq fJ^ + 6(1 - a)«;(/ - a^A {I - a) + [as, r]p, (128) 
D_g(/ -a) = -qas^I - cr) + 2 fJ^ - uO) (/ - a)pas - q[as, r], (129) 



n _ _i 



I?_p(/ + a) = ^^3^(1 + a) - 2(/ + a) ( W- + (/ + a)-^u;(/ + a)e 1 ^^3 - ^[^3, r]. (130) 

Equations (I109p . (IllOp for the derivatives of r together with a consequence of (11181) (since 
dkR = dkp) [23], 

dk^ij = dk{Rij{^,^j)) = idkRij){^i,^j) + dxRij{iijij)5ik + dyRij{^i,^j)6jk = 

= -TikTkj + dxRiji^i, QSik + dyRij{^i, QSjk, (131) 

give (we will sometimes write T>^ for (cr30")I?_) 

P+r = r^ + ry-r^ = -p{I + a)Qq - qQ{I - (r)p - [a^, r], (132) 

V^r = ar^ + Vya — rar = —p{I + cr)9g + g9(/ — a)p — [(,, r]. (133) 

Besides, differentiating the definitions (11121) . (I113p and using flUTl) . (jHSD and flU^ . we get the 
matrix analogs of universal equations [2T] for scalar u and w: 
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V^u = —qq, V^u = —qa^^q, V^w = —pp, V^w = —pa^p. 

Thus, we have obtained a system of fourteen first-order matrix PDE for joint largest 
eigenvalue distribution of two coupled Gaussian matrices. Our system is different from that 
in [23] since we employ the matrices u and w. As we will see, this system is again the most 
convenient for comparison with r-function approach of Adler and van Moerbeke [3l H] like it 
was the case for one-matrix ensembles [T71 ITS] . 

We remark that we are interested in the matrix r while equations fll32p . fll33p involve 
only terms where the other (auxiliary) variables enter multiplied by the degenerate constant 
matrix 0. Therefore, since B^ = 20 and 



{I-a){I + a) = I-a'= l-'-—^\.I=-—-.I, (134) 



[i + cyj (i + c)2 

it is convenient to introduce new matrix variables, which in fact amounts to some reduction 
of the total number of scalar variables since the new matrices, being proportional to O, also 
have determinant zero: 

i^,e, l^ei, p.(i^fML±^e, p^e ^-f^;^\ (135) 

1 — 0"^ \ — a'^ 



1e-e.„e. H/ . w„ = ^fe + e(ii:^^>^^;ii±^a (lae) 

Here and further on, in a slight abuse of notation, we write a"^ for the scalar (1 — c)^/(l + c)^, 
so 1 — cr^ stands for 4c/ (1 + c)^. To express everything in terms of these new variables, we 
multiply equations (I123p . (I124p . (I127p . (I128p by the matrix 6 on the right and equations 
dnSD, dnSD, dnnD, (HSOD by the matrix e on the left. This brings our TW-type system to 
the form given in theorem |U thus finishing its proof. 

It is easy to verify that the system has two matrix first integrals, similar to the first integral 
of [21] for Gaussian single matrices: 

Pq = nQ- WU, (137) 

^p = ne- UW. (138) 

6 Coupled analogs of one-matrix relations among ma- 
trix kernel related variables and r-functions 

There are direct analogs of one-matrix universal relations [IB] here. 
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Lemma 2. Analogs of rank-1 projection operator kernels read: 

Kn+iCL - eiKn = eLQ(Pn{x)ipn{y)x{y), (139) 

euKn+i - KnCu = Qeu^n{x)Vn{y)x{y)- (140) 

Proof. The diagonal entries of the matrix kernel K are the same as in scalar (1-matrix) case. 
Consider the anti-diagonal entries of the kernels for matrices of consecutive sizes n and n + 1. 
Since 

n— 1 CXI 

{Kn)2i{x,y) = Y^d'-^^k{x)^k{y)x{y): {Kn)i2{x,y) = -^c^->fc(x)(^fe(t/)x(y), 

A;=0 k=n 

we find 

-(^n+l)21 - {Kn)2l{x,y) = ipnix)ipniy)xiy), c{Kn+l)l2 - iKn)l2ix,y) = V?„(x)v?„(l/)x(Z/), 

c 
which means, that we have matrix equation 

-euKn+iix, y)eL - K^ix, y) = Qipn{x)ipn{y)xj{.y)- 
c 

By the identity eiCu = cuCl = cl, it is equivalent to the statement of the lemma. 

D 

Lemma 3. There are the corresponding relations for the resolvent kernels: 

Rn+lCL - ClRu = Pn+l{x)eLQQn{y), (141) 

euRn+i - Rneu = Qn{x)QeuPn+i{y)- (142) 

Proof. Writing Kn+iCi - eiKn = ciil - Kn) - (/ - Kn+i)eL and multiplying f ll39p by 
I + Rn = {I — Kn)~^ on the right and by J + Rn+i = (/ — Kn+i)^^ on the left gives equation 
flTiT]) . Equation flT42|l follows from firiOjl similarly. 

D 

Also by definition of matrices Qn, Pn, Qn and P„ we have recursion relations: 

ifin = {I- Kn)Qn = (/ - Kn+l)Pn+l = Qn{I - K^) = P„+l(/ - K^+l) . (143) 

They in turn lead to relations involving inner product matrices u„ = {ipn, {I — Kn)^n) and 
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Lemma 4. 

6(7 -Pn = Qn-l(^U + Q„-10e(7W„, -Pn+lCi = ClQ^ + Wn+ieL^Qn, 
eiQn = Pn+ieL - Pn+ieL&Un, Qn-lGu = (^uPfi " Un-lQcuPn- 

Proof. Let us prove e.g. the first of the above formulas: we have 

'^n^U = eu{I - Kn)Qn = eu{I - Kn+l)Pn+l = {I - Kn)euPn+l " {tuKn+l - -^net/)-P„+l = 

= (J - Kn)euPn+l - (fn&euWn+l, 

and acting on both sides by (/ — Kn)^^ from the left gives 

i.e. the sought formula if we shift n — )■ n — 1. The other three formulas are obtained quite 
similarly. 

D 

Introduce matrix w„ = ejjWneL/c, then we have 

Corollary 9. 

I + Qijn+l = {I - QUn)-^ (144) 

and 

I + ^n+iQ = {I - UnQ)-' . (145) 

Proof. One can rewrite recursion relations from previous lemma as follows: 

— = Qn{I + QWn+l) = Qn{I " 0M„)"\ 

C 

^uPn+lGL 



C 

thus we get the statement. 



(/ + w„+i6)(5„ = (/-M„e) ^Qr 



U 



28 



It follows also for the matrix resolvent kernels that 

Rn+ieL - eiRn = eLQn{x){I - QUnY^QQniy), 
euRn+1 - RnCu = Qn{x)Q{I - UnQY^Qn{y)eu ■ 

Taking x = ^j, y = ^j in the last formulas gives, respectively, 

euTn+l - TnCu = ?„©(/ " U^QY^ Qneu ■ 

Now, after recalling again eiCu = ejjGL = c/, it follows from either of these two formulas that 



V+\n 


' n 


= Tr(r„+i - 


-Tn) 


which can 


be integrated to 


give 








'^n+l/'^n+l 


- H 



^ , det(/ - Qun) = det(/ - M„e) = 1 - Tr(eM„) 

Using the recursion relations to get the second, similar formula, we thus prove the main result 
of this section - theorem |5l since we also get the formulas for the matrices U and W defined 
before it in section [3j 

Trf/ = J-TrU = \pM - 2Tr(0M„) = v^(l - Tr(0M„)) = v^ '^n+Jju+i ^ 



Tr# = J-TxW = V2^ + 2TT{ewn) = v^(l + Tr(0^„)) = v^ 



In fact, one has a remarkable result: 



Lemma [T]. 

Tr{UW) = TiU ■ TvW (146) 



29 



Proof. Since Trr = V^lnr^, Tr(cr3r) = V_\nT^, Tr(o"r) = V^lnr^, we get from the equa- 
tions ([121 and (1^5]) : 

We use the first integrals ( I137p and f ll38p to find 

Tr(gp) = Tr(pg) = Tr(ne - f/#) = 2n - Tt{UW), 
so the previous equation can be written as 

(Vl - Vl) Xnri = 2(1 - a') QTr(?7#) - n}j . 

We compare this last equation with the "boundary- To da" equation from section 3, which can 
be written as 

{Vl - Vl) Inri = 2(1 - a2)(2(l - ^)UnW^, - n), 
which immediately leads to 

Tr(f/#) = 4(1 - C^)UnWn. 

Therefore, by the previous theorem, see also section 2 for the normalization of the 2-matrix 
integrals r„ over the whole domain, 

Trf/ • TrW^ = 2nUnWn = 4(1 - c^)UnWn. 

'Tn+l'Tn-l 

n 

7 Transformation of the TW-type system and PDE analogs 
of Painleve IV 

We consider and solve the system obtained in theorem O 

Lemma 5. The system of theorem\^ has three full matrix first integrals: 

[X+,G] = {X_,$}, (147) 

[X+,<I>] = {X_,G}, (148) 

[$, G] = {X_, 3X1 - 8nX+ + X!} - [X+, [X+, X_]], (149) 
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Proof. The first two formulas are just tlie integrals fll37p and (11381) rewritten in new variables, 
and the third one can be easily derived from the above system of PDE as consistency condition, 
if using the first two. Consider e.g. 



[P+X+, G] + [X+, V+G] = [$, G] - [[e + A, X_], G] + [X+, {^ + i, $}] + [X+, [X+, X_]], 



I.e. 



[$, G] = P+[X+, G] - [X+, [X+, X_]] + [[e + i, X_], G] - [X+, {e + i, $}] = 

= P+{X_, $} - [X+, [X+, X_]] + [[e + i, X_], G] - [X+, {e + i, $}] = 

= {X_, 3X2 _ g^^^ + X' } - [X+, [X+, X_]] + 

+ [[e + i, X„], G] - [X+, {e + A, $}] - {[e + i, X+], $} + {X_, (e + i, G}}, 
which gives (11491) . since the last four terms on the last line cancel out due to f ll48p . 



D 



There are additional diagonal matrix first integrals, which can be obtained by integrating 
the diagonal parts of equations ( l38l) and, using expression for $ from eq. (1371) . also ( HOl) . or 
diagonal parts of equations fH2l) and (jH]), using expression for anti-diagonal part of G from 
eq. f H2|) on the right-hand side of (jS]). The pairs of equations obtained in these two ways 
become identical after integration. Thus we get, respectively. 

Lemma 6. The system has two diagonal matrix first integrals, 

(X_), = -AV, (150) 

Gd = 4r, - 2^+Vrd - 2e_P_r, + {A, (X+)4 - a-siA, XJ. (151) 

Splitting into diagonal/anti-diagonal parts appears to be convenient for the rest of the 
system also. Equations for the derivatives of r in diagonal/anti-diagonal splitting, if we 
denote 

are 
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= -Xd + a(X_)d, 2P_rrf = -aX^ + (X.)^, 



2D+rd = -Xd + a(X_)d, 2P_rrf = -aX^ + (X.)^, (152) 

asV+A = -{X+)a - ^+A, aV^A = X, - ^-A. (153) 

The equation (I150p then entails important simple relations 

X, = -2V+rd - a^A^ = -2a;V^u - ^'. (154) 

Further splitting of the diagonal parts into scalar trace, we denote TrM by Mt, and tr(o"3M), 
denoted by M3, parts will also be used. As it is clear from the eq. f l38|) . TrX_ = 0. As follows 
from the connection of the matrix r with Inr;^, 

(the second equality above can be seen also from each of the equations (11521) ). Also, a 
consequence of (11541) is the formula [26] simply relating anti-diagonal elements of r with its 
diagonal elements. 

One can see that, besides equations (I132p . (I133p . some combinations of the other equations 
in the system are also linear in X_|_,X_, $ and G. Namely, adding f l57|) and (crqa)- (H21) gives 

v^x, + Mv^x^ ^ * - ,JJ^ - ,,M±1 , Ml. 

while adding ( 138|) and (o"3(t)-(J4T]) results in 

It follows from the linear equations that 

^t = V+Xt = P_X3 - 2V+A\ (155) 

$3 = V_Xt = V+X^ - 2a'^V_A^, (156) 

recall that a^ = (1 — c)^/(l + c)^. Recall the quantities A^, A_, X3, Gt, G3 defined in the 
corollary [3, in section [3l 

The splitting of eq. (I15ip into the scalar parts reads: 

Gt = An - 2i+V+n + ^^3 + A+ = Ht + A+, (157) 

G3 = 4r3 - 2e-P-r3 + ^+^3 + A_ = H:, + A_. (158) 

In fact, one can show that 
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Theorem 10. The system can be reduced to the following five independent (scalar) ^th-order 
PDE mhir!: 



where 



X^V+V^X^ - P+XaD-Xs + RtR-^ - (X3 + i+^-)Xl = 0, (159) 

Rt = Art - 2^+I^+n, i?3 = 4r3 - 2i_V_r^, 

V+V.Xt = i-Gt + e+Gs + X3(3X, - 8n), (160) 

2FX^VlXt = 2F{V+XsV+Xt - RtG, + ^+X,Gt) + Xs{V+X^ - G?), (161) 

2FX3V^_Xt = 2F{V_X^V.Xt - R^Gt + ^-X.G^) + X^iV^X^ - Gl) (162) 



and 



'''^'^n-a^^' = -2a2l?+P_A2 = -{Qa'XsA' - ^_A+ - a'^^A^). (163) 

For the proof, see Appendix B and formulas 0185p and 0186p below. 

It is rather tricky to integrate these equations directly. There is, however, a simpler way. 
Additional integrals can be most readily seen from the original defining variables: 

det(g) = det(g) = det(p) = det(p) = det(f/) = det(iy) = 0, 
so 

det(X+ ± X_) = det($ ± G) = 0. (164) 

There are also remarkable identities here: 



Lemma 7. 



Tr ({UW)'') = Tr ({WU)'') . (165) 



Proof. It follows by induction from the case k = 1 proved in the previous section by compar- 
ison with the results from Toda lattice. 

D 

Then we retrieve 
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Lemma 8. It follows from eq. (TU^ that (besides eq. li20S\) in Appendix B which arises this 
way also) 

4((X+)^ + Xl) = J^Xf- Xl - Aa\A^)\ (166) 

ml + Gl) = $? - $1 + G? - Gl (167) 

2{<l>a,Ga} = Gt<^t-Gs<^s, (168) 

and from U65\) one can get once again both eqs. li20S\) and U66\) . but also another first 
integral: 

We will need three more auxiliary formulas, which are direct consequences of previously found 
first integrals fl200D . fl20TD in Appendix B: 

XlGl = Gl{X^)l - <I>?X,2 - $,^3^, (170) 

X|$^ = ^liX^)l - G\Xl - $3^,^, (171) 

X|{$„ Ga} = 2($3G'3(X+)^ - <^tGtX'J - ($^$3 + GtG,)Ga. (172) 

First two of the last equations are obtained by taking square of eqs. (12001) and (I20ip . respec- 
tively, and the third is their anti-commutator. Plugging (I166P into eq. (12081) of Appendix B 
simplifies it, giving 

$^$3 - GtG, = 2X3(X2 - 4nX,), (173) 

Using these relations together with lemma 8, one can verify that all the higher-order equa- 
tions are satisfied, see Appendix C. 

One can easily eliminate all auxiliary variables but two: A^, A^, and get 

Theorem 11. The joint largest eigenvalue distribution for two coupled Gaussian matrices 
satisfies the system of PDF: 



$^$3 - GtGs - 2x,X3F = 0, mm 



A+A_ = Aa\V+A^V_A^ - 2X3(^^)2), (174) 
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($t)' + ($3)' - G? -Gl- 2F{2Xl + J) = 0, (175) 

A\ + a'^A^_=Aa'^{Dl + D't + A^J), (176) 

F{Al - Aa\Dl - XlA^)) = 2o''A^{^\ - G\), (177) 

F(A^ - 4(a'D2 - XlA^)) = 2A^{^1 - Gl), (178) 

F{D+A^ - D_A+) = A^^tGs - $3^0, (179) 

with notations introduced: 

F = Xt-An, J = X^-Xl-Aa\Ay, D+=V+A\ D^=V_A^. 

For the convenience of the reader we write out again the definitions of the variables entering 
the system above: 

Gt = Ht + A+, G3 = Hs + A., Ht = Art-2^+V+rt + ^-X3, H3 = Ar3-2^^V^r3 + ^+Xs, 
$i = V+Xt = P-X3 - 2P+A2, $3 = V^Xt = P+X3 - 2a^V^A^ a^ = (1 - c)V(l + c)^ 



X3 = -2P+r3 = -2V^rt = -2VV_ In r, 



Proof. We already have equations (I173P and (I174p (for the last see lemma 11 in Appendix 
B), equation (11761) is the result of taking the combination of equations (I210p + cr^ (12 lip from 
Appendix B and using (I166p to eliminate {X^)"^^ + X^. After expressing G^ and $^ from 
equations (I167p and (I169p : 

AGl = ^l-Gl- 2X^F, (180) 

4$2 = -$2 + G^ + 2X^F, (181) 

and plugging them into (11701) and (I17ip . respectively, the last become: 

Gl{A{X+)l + Xl) - ^l{AXl + Xl) = A^,G,Ca - 2XlXlF, ^M 
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$2(4(X+)2 + Xl) - G'.iAXl + X|) = A<^,GtCa + 2XlX^F. (171]) 

Applying ( IMI]) . (12UB and (IT75]) transforms equation (IT72D into 

$3G'3(4(X+)^ + X|) - $tG'i(4X2 + X|) = 4(<l>t<l>3 - X^XtiXt - An))Ca. mM 

Then we make combinations, $f ■ f ll72p — Gt ■ ( I170p . which gives, after using fll73p . 

2X^XtFGsmX+)l + X|) = AXsXtF^tCa + 2XiX^FGt, 
and -$3 ■ (IT72|) + G3 ■ (IT7T]) . giving, after (IT73|) is apphed, 

2X3XiFGi(4X„^ + X|) = 2XlX^FG3 - AX^XtF^sCa- 
They are obviously simplified into 

G'3(4(X+)2 + Xi) = 2<^tGa + X.XtGt, (182) 

Gt(4X,2 + X|) = X,XtG3 - 2<^sCa. (183) 

Their combination $3 ■ (I182p — $< ■ (I183p . compared with (I172p . means that 

4<l>,$3a + X3X,($3Gi - <^tG,) = 4($,$3 - X3X,(X, - 4n))Ca, 
i.e. 

$,G3 - ^sCt = 4Fa. (184) 

Putting expression for Ga from equation (12121) of Appendix B in (11840 proves (I179p . Mul- 
tiplying (I182P and (I183P by 2F and replacing 4FC„ in them by the left-hand side of (11841) 
yields 

2F(4(X+)^ + X|) = $,2 - G? = V+X^ - Gl (185) 

2F(4X,^ + X|) = ^l-Gl = V^X^ - Gl (186) 

Adding (11851) and (11861) and using (11661) gives eq. (11751) . while using them to eliminate (X_|_)^ 
from eq. (I^TUD and Xf from (I^TTD leads, respectively, to eqs. (ITTTp and (IT7g]) . D 

The system of equations obtained is still redundant since equations ( ]175p . (11760 . ( 1177P and 
(I178p are in fact linearly dependent. Therefore The system of PDE in the previous theorem 
is equivalent to the one in the main corollary [71 
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8 Comparison with Toda lattice approach and new re- 
lations among different variables 

Let us compare with the results of Toda lattice approach in section HI Recall that there we 
used functions 

F = UW = UnWn{l~c^), G+ = WV+U -UV+W, G^ = WV_U -UV^W. 

Certain combination of equations f ll60p of theorem [TD] from the previous section and (11731) of 
the main corollary [7] is in fact exactly equivalent to the equation ( 189|) from theorem [3] derived 
in section m The combination is 2-F- fll60p — f ll73p as can be guessed by comparing its senior 
derivative terms of 4th and 3rd order with that of (!89|) . It can be written as 

2FV+V_Xt - V+XtV^Xt + GtG^ - 2F{^+G^ + ^GJ - 4X3^^ = 0, 

which clearly has the same form as (1891) . Comparing the terms in the two equations, we find 
simple correspondences: 

4F = ^^^^-^'f + 2- = -1/2(X, - in) = -1/2F, 

G-(- ^ Gt/o, G = G-3/8, 

A+ = 4((1 - c^)d,n - i+<y^A^), A^ = -4((1 - c')d,r, + ^_A'), 



and so 



asK V+ra] = -2cdcrt, ds^, P_r,] = 2cdcr3. (09]) 



These last two new important relations show that the main additional variables appearing 
in the TW approach to coupled (Gaussian) matrices - the commutators on the left-hand 
side - directly correspond to the main new objects appearing in the ASvM approach to the 
same problem - the derivatives w.r.t. the coupling c (or w.r.t. time t, recall c = e~*) of the 
log-derivatives of the joint largest eigenvalue probability. 

The Adler-van Moerbeke equation [3], formula (182|) of section HI in the current variables 
reads: 

,« G'i ,„ Gf — 

F F 
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9 The smallest complete set of PDE for two Gaussian 
coupled matrices 

One can in fact eliminate the remaining auxiliary variables - A^ and A_, which turned out 
to be directly related to the derivatives of In t^ with respect to the coupling parameter c, and 
thus explicitly obtain PDE in terms of spectral endpoints only. To this end, let 



Ft = F + 2A^, F3 = F + 2a^A^, A = FtH^ - F3HI J = X^ - X^ - Aa^{A^ 



J+ = 2a\2Dl + A^J), J^ = 2{2a^D^_ + A^J), Ja = Aa\D+D^ - 2X3{A^f), 

P, = ^t^3-HtH3-2X3XtF-JA, Pt = ^l-Hl-{2Xl+J)F-J+, P3 = ^l-Hl-{2Xl+J)F-J, 

P+ = FtPt + F3P3, Pa = HlPt + H^,P; - 2HtH3P.. 

Recall that 

F = Xt- An, D+ = V+A^, D_ = V_A^, 
and introduce 

St = A^^t - FD+, Ss = A2$3 - FD_, Ja = Hs^t - Ht^-,. 

Then, expressing y4+, A^ from the equations in the main corollary [TJ which are linear in 
them and their squares, and putting into the other ones, one gets four final equations: 

{HtP^ - 2F:,HsP:,){2FtHtP, - H3P+) = 4A'Ja, (187) 

(HtP+ - 2F3if3P.)' = 4A(AJ+ - 2a^A^Pa), (188) 

{2FtHtP. - H^P+f = 4A(AJ_ + 2A^Pa), (189) 

S^iHtP^ - 2F,H,P,) - St{2FtHtP^ - H,P+) = 2A^AJ,, (190) 

as one should since there are four independent senior derivatives involved here: 

38 



$i = V^Xt = V_X3-2V+A'^, $3 = V_Xt = V+X3-2a'^V_A'^, D+ = V+A^, D_ = V_A^. 
From the first three one derives: 

Pl = P^ + 2{H^,-a^Hl)^, (191) 

Pi = A{F,F,Pa + A(F,J+ - F,J_)) + %{FlHl - a'pfH^)^, (192) 

A^Pl = 2APaiDl - a^Dl) + A^/^, (193) 



where 



Pa = HlJ+ + H'^J_ + 8a'^HtHi{D+D_ - 2X3(A 



2^2^, 



la = J{2{Dl + a^D^ ) + ^2 J) ^ i6a2A2X3(D+D_ - X3(A2)2). 
We repeat the involved definitions once more for convenience: 

Xt = -2V+rt - 2a^/^ = -2V^r^ - 2A^ , X3 = -2V+r3 = -2V^rt, 
n = V+hiTi, r3 = P_lnr„^ A^ 



l-a2 ' 



Ht = An - 2^+V+n + ^_Xs, Hs = 4r3 - 2e-I^-r3 + ^+^3. 

The combinations Px and P+ turn into Painleve IV equations in the one-matrix hmit, while 
Pa, Pa, la as well as A go to zero then. 
Since 

A = FtHf - FsH^ ~ 2(1 - a'^)A'^{8Gf, 
and, when ^2 -^ 00, 8G = 4(rt — (.id^j^rt) = 4(r3 — ^19^1^3), we can estimate 

Pa = HlPt + H^P; - 2HtHsPx ~ {8Gy{CiA^ + ^2%^^), 
where the quantities Ci and C2 remain finite in the next limit, so at fixed c (or a^), as ^2 -^ 00, 

A^Pa 

— r^-^0. 

A 

Thus, equations 01911) and (11921) also tend to Painleve IV in the one-matrix limit, while the 
other two - (J193I) and (I19UI) - become trivial. 
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10 Conclusions 

The joint probability for the (largest) eigenvalues of two coupled Gaussian matrices with 
unitary invariant probability density satisfies a number of nonlinear integrable PDE some of 
which are coupled analogs of Painleve IV equation for one-matrix GUE. The corresponding 
equations for the scaling limits of Airy process can be obtained from the system in the main 
corollary[7]of section[3l but their consideration as well as due generalizations to several spectral 
endpoints (partly treated in section H] though) and several coupled matrices are delegated to 
a forthcoming work. 

The matrix kernel approach appears in a sense superior to the one based on Hirota bilinear 
identities and Virasoro constraints. The first allows to obtain at once all PDE satisfied by 
the joint gap probabilities while the last gives various subsets of the whole set of such PDE, 
depending on which nonlinear integrable equations in "times" have been taken as the starting 
point. This situation is quite similar to what we recently found for the one-matrix case with 
several spectral gaps [19]. 

For the Gaussian coupled ensemble simple relations among variables of different approaches 
are found here, and there is hope, supported by our analysis of single-matrix UE [TS], that 
they again can be extended to other coupled RM. For some of them this is already shown 
here in section [61 Comparison of biorthogonal function structures considered in [^ and the 
matrix kernel approach of [23] is an interesting future direction. Biorthogonal case analogs of 
3-term relations for functions in the resolvent kernel [IB] , which may combine finite difference 
recurrence and differentiation formulas, need to be found. They would present the coupled 
case from the (bi) orthogonal functions point of view, appearing the best for the description 
of one-matrix ensembles. 

A Dyson Brownian motion (BM) model (or Dyson process) for all orthogonal-polynomial 
matrix ensembles can be constructed. It gives a class of coupled ensembles whose joint prob- 
ability density satisfies a diffusion Fokker-Planck (FP) equation. The BM model can be 
mapped onto a quantum many-body problem and its transition density conveniently repre- 
sented in terms of Green function for the quantum mechanical model, see e.g. [8l [13]. Then, 
for unitary ensembles possessing the BM representation, one can get the matrix kernel of the 
form obtained in [23] for Hermite and Laguerre cases, i.e. express it in terms of the corre- 
sponding eigenfunctions and eigenvalues of the effective Hamiltonian obtained from stationary 
FP operator by the above mapping. Then it is possible to give a TW-type derivation of PDE 
for arbitrary BM ensembles, if the differentiation formulas [21j are available, a situation just 
like in the one- matrix case. 

Current work also paves the way for obtaining new integrable PDE satisfied by gap prob- 
abilities of other ensembles with similar matrix kernel such as e.g. Pfafiian ensembles. 
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Appendix A 

Proof of theorem [2] is given here. 

For further analysis we need more commutation relations: 

Ao- = -{Ao + 1) - 2A, Ao- = -{Ao + l)- 2A. 
c c c c 

The system of five equations we are considering apparently contains more unknowns than 

equations but we can enhance it. To this end, add up eq. ( !66|) multiplied by AW and eq. ( !67l) 

multiphed by AU, which leads to the relation: 

2{AWAoU - AUAoW) = -2A^TAF - A ( ^^^IrT^ ) • (1^^) 

The "dual" copy of it is the outcome of adding up eq. fl68|) multiplied by AW and eq. fl69l) 
multiplied by AU: 

~ - ~ ~ ~ f (AF)^ - G^\ 
2{AWAoU - AUAoW) = -2A^TAF - A ^ ' . (195) 

The expression on the left-hand side of fll94p can be written in two different ways: 



AWAqU - AUAoW = AGo - {WAAqU - UAAqW) = -AqG + {WAqAU - UAqAW). 
Using the commutation relation 

AqA = A{Ao - 1), 
one gets 

2{WAAoU - UAAoW) = {Ao + 1)G + AGo, 
and so 
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2{AWAoU - AUAoW) = AGq - (A + 1)^. 
The "dual" of this also obviously holds: 

2{AWAoU - AUAoW) = AGo - (A + 1)^. 

Plugging the last expressions into the left-hand sides of eqs. fll94p and fll95p . respectively, one 
gets the needed additional equations 

AGo = (A + 1)^ - 2A'TAF - A ( ^'^^\~^ ) ■ (196) 

AGo = (A + l)G - 2A^TAF - A I ^'^^\~ ^ ] ■ (197) 

Using the found expressions for G, G in terms of T and commutation relations, we can express 

AT AT 

(A + i)G = -(A + 1)(A - 1)^ = HAl - 1)^ = 



= -i/2(A(^/c(A+i)-2A-^/c)T = -i/2((Ac(A+i)-2A(A+i)-2^(A-i)-^/c)T 

= -l/2{A/c{Al + 2 A) - 4^A)r, 
i.e. 

(A + 1)^ = 2AAqT - ^A{Al + 2Ao)T. (198) 

Completely similarly, 

(A + 1)^ = 2AAoT - ^A{Al + 2Ao)T. (199) 

The last expressions together with ( !74l) mean that the eqs. (11961) and (I197P can be written as 

A (go + ^^^\~ ^' - 2At] = -A (\-y^Tf + ^yl + 2A)t] , ([MD 



^ I Go + ^^^\p ^' - 2AT j = -A {^^{A^Tf + ^(^0 + 2A)T ) . (IMD 
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Now we eliminate Gq and Go from them with the help of eqs. (1771) and 0781) and obtain two 
higher-order, in fact 5-th order in T, PDE where everything can be expressed in terms of T 
only, the ones in the theorem [2l 

Appendix B 

Let us write out the diagonal and anti-diagonal parts of the matrix equations fll47p . (11481) 
and (I149P separately. Equation fll47p gives, using fllSOp . 

which splits into two scalar equations: 

[(X+)„ Ga] = -$tAV, {Xa, $a} = ^S^'l^S^). 

Its anti-diagonal part reads: 

X,Ga = G3(X+), + <^tasXa. (200) 

For the eq. (I148p we get in the same way: 

[{X+)a, $„] = -2AVG, + {Xa, Ga}, 

which splits into 

[{X+)a, $J = -GtA^a, {Xa, Ga} = G,A\ascj), 
and the anti-diagonal part, 

Xs<^a = $3(X+), + GtasXa. (201) 

Combining the above equations implies some new ones, e.g. 

{iX+)a,Xa} = iasa)X3A\ (202) 

^3[$a, Ga] = {GtGs - $^$3)^^, (203) 

X.i^sGa - Gs^a) = ($t$3 " GtG,)a;Xa. (204) 

Taking the diagonal part of equation (I149J) (using (11501) again) gives two scalar equations, one 
of which is the above (12021) again, while the other is 
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3X^ + X^ , n^v \2 , nv2 \ a2. 



l^a, Ga] = - ( ^ '- - 8nXt + 2(X+)^ + 2XM A'a. (205) 

The anti-diagonal part of f ll49p is 

$3^3Ga-G3a3$a = (3trX2-8nX,+2X!)X,-X3a3(X3a3X,+2AV(X+)J+2[(X+),,X,](X+),, 
i.e. 

<^3Ga-G,<^, = I '^ ' - 8nXt + 6(X+)2 + 2X!j a3X„+2(a3[(X+)„X,]-X3AV)(X+),. 

(206) 
One has an identity for 2x2 anti-diagonal matrices: 

2(X+)2X„ = ({(X+)„X4 - [(X+)„Xj)(X+)„, (207) 

which, together with (I202p means 

2(X+)2a3X, - (X3AV - a3[(X+)„Xj)(X+), = 0. 
The last relation simplifies equation ( 1206p . it becomes 

$3G'„ - G'3<l>. = i '^ ' - 8nXt + 2(X+)2 + 2X!j agX,. mB 

Comparing this with eq. (I204p one gets a scalar relation 

$,$3 - GtGs = X3 '^ ' - 8nX, + 2(X+)2 + 2X! , (208) 



which equally well follows from eqs. (I203P and (I205p . There is actually a lot of redundancy 
here and some relations arise in several different ways. 
To summarize, we proved 

Lemma 9. The full matrix first integrals reduce to four independent relations - two scalar 
and two anti- diagonal matrix: 



{{X+)a,Xa} = ia3a)XsA', (ESD 

X3G, = G3(X+), + <l>,a3X„ (HP 
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Xs^a = ^3{X+)a + GtasXa, ([MI) 

$<$3 - GtG^ = X3 * ' - 8nXt + 2{X+)l + 2X! , mM 



and all the other above relations follow by taking commutators or anti- commutators involving 
mm and/ or mi\) . 

Next we use the identity (of which eq. (12071) is an instance) 

Lemma 10. For any anti-diagonal 2x2 matrices Aa and Ba'- 

2Al- B,= {{Aa,Ba} ~[Aa,Ba])Aa 

and its scalar consequence 

AAlBl = {Aa,BaY-[Aa,Baf. 

In particular, matrices (X+)a and Xa = {X^)a can be expressed as 

2A\X^)a = i{A,X^}-[A,X^])A. 
Besides, we have the facts: 

[A, (X+)J = P+AV3, [A,Xa] = -V_A'a, {(X+)„XJ = {a,a)X,A', (209) 



first two are consequences of the eqs. (11531) and the third is the eq. (I202p . We apply all this 
to derive 

Lemma 11. 

16a'^A^{X+)l = A\- Aa'^{V+A^f, (210) 

l&A^Xl = A^_- Ao'^{V_A^f, (211) 

AA^Ca = 4AV3[(X+)„XJ = -A+V_A' + A_V^A\ (212) 

A+A^ = A(x^{V+A^V^A^ - 2X3(^2)2). ( [77^ 



The last equations should be considered as four more scalar first integrals. 
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Let us substitute the expressions fl200p and (12011) into the right-hand sides of the equations 
of the system. The anti-diagonal parts of the equations, besides f ll53p . become: 

X3l?+(X+), = <I>3(X+), + {Gt - ^-X^)a^Xa - 2X3AVi, (213) 

a;V+Xa = X,A-i^{X+)a, (214) 

X^V+^a = (e+Ga + X3(3Xi - 8n))(X+), + e+^to^s^a + ^sG^i, (215) 

X^V+Ga = i+^z{X+)a + 2X3AV(X+)„ + {i+Gt + XDa^Xa + X^^.A, (216) 

V_{X+)a = -Xsa-iA - i+a^Xa, (217) 

X^a^V^Xa = (G3 - e+X3)(X+), + ^t^y^Xa + 2X3^2^, (218) 

X3P_$„ = (e-G3 + X|)(X+), + (^<l>i + 2X3AV)(T3X, - X3<l>3a3A. (219) 

X^V.Ga = e-$3(X+), + (e-Gt + X3(3X, - 8n))a3X, - X3G3a3A, (220) 

Applying the above anti-diagonal identities simplifies eqs. (I213P and (I218p . they turn into 

X^V+{X+)a = V+X^{X+)a + {Gt - e-X3 - A+)a^Xa = V+X^{X+)a + RtCX^Xa, ^ 

X3(T3l?-(X_)„ = (G3 - e+^3 - ^-)(^+)a + ^-X3(T3X, = R^{X+)a + I?+X3a3X„ dm 

where 

Rt = An - 2i+V+n, R, = 4r3 - 2^P_r3. (221) 

As follows from these definitions, 

V_Rt = e+I^+Xg - 2X3, P+i?3 = e-^-^3 - 2X3. (222) 

We now match the mixed second derivatives T'+P_(X+)a etc. found from the corresponding 
pairs of the above anti-diagonal equations. It turns out that taking both pairs of equations 
), dnZD and (EH, dnSD, as well as plugging ^M^ and ([MD into the left-hand sides of 
and (12191) . respectively, leads to the same higher-order scalar consistency equation: 
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Lemma 12. 

X^V^V_X^ - V^XsV^X, + RtRs - (X3 + e+e-)X| = 0. mB 

This is in fact a 4-th order PDE for liar^, since all the variables in Iil59\) have been already 
expressed in terms oflnr^ before. 

However, we are going to obtain 3-rd order equations in Inr;^ and show that (11591) is a 
consequence of them. 

Consider next the eqs. (12141) and fl217p . Applying (11531) . one sees that both of them then 
reduce to the same anti-diagonal consistency equation: 

P+P_A=(X3 + e+e-)A (223) 

Using it, compute 

V+V^A^ = {A, V+V^A} + {V+A, V^A} = 2(X3 + ^+C-)A^ + {V+A, V^A}, 
and, applying (I153p again and 



a^{V+A,V_A} = -l/2^+a^A_ - l/2^-A+ + (X3 - 2^+e_)aMl 
Equation (I163p of theorem [TU] immediately follows from the last two formulas. 

Lemma 13. The remaining anti-diagonal equations Ii215\) and Ii220\) . if we substitute Ii201\) 
and Ii200\) . respectively, into their left-hand sides, result in three independent scalar equations 
(three rather than four because T'+$3 = T>_^t = V^V^Xf): 



I?+<l>3 = P_$i = ^-Gt + ^+Gs + Xs{3X, - 8n) 



XsV+Gt = V+X,Gt - Rt<^s + e+X3$t, (224) 

X3P_G'3 = V^X.G, - Rs^t + e-X3$3. (225) 

The last two of them, however, have been in fact already integrated when Gd was obtained: 
they are trivially satisfied if Gd is plugged into their left-hand sides. 

The other 4-th order equations can be produced in more than one way. We write out just 
one here. 

Lemma 14. The following higher-order (4-th order in lnr;^j equations hold: 

X^V+^t = P+X3$i - RtG^ + i+X^Gt + X;{Xl + 4(X+)^), (pS!) 

X3l?-<l>3 = I?-X3$3 - RzGt + ^-^3^3 + X^{Xl + AXl). (USi 
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Proof. Consider the remaining diagonal equations of the original system, namely the ones for 
the derivatives of $d: 

!?+$, = 3(Xj + {X+)l) - 8nXd + X^_ + 2^Gd + {i, GJ, (226) 

V^<^, = 2a,{Xl + {X+)l-AnX, + Xl) + a,{Xl-{X^)l)-a,{{X^)l~Xl) + 2a,iGa+[A,^a]. 

(227) 
Substituting ( IMJ]) into (I22B]) and (12UB into (12^ . then separating the Tr(-) and Tr((T3-) parts, 
we get equation (11601) again plus two more scalar equations: 



X^V+^t = X, { —^ i - 8nXt + 2{X+)i + 2Xt\ + X^{X', + 4(X+)f ) + 

+ i+X^Gt + (e-Xs + A+)G; + 2a2l)_A=^$,, (228) 

8nXi + 2{X+)l + 2X!') + X^^^XJ + 4X2) + 



X3P_$3 = X3 



3X,2 + X| 



+ (e+X3 + A^)Gt + ^-^3^3 + 2P+A2$3. (229) 

Now we recall relations (I155p . (11561) and (12081) and see that these last equations reduce to 
equations (I16ip and (I162p . respectively. 

D 

This, together with formulas (I185p and (I186p applied in eqs. (I16ip and (11620 . ends the 
proof of theorem [TOl 

Appendix C 

In Appendix B we proved theorem [TO] - the five 4th-order equations for Inr;;^. 
After using formulas 

2F{A{X^)l + Xl) = <^l-Gl dUSD 



2F{AXl + X|) = <^l- Gl ^ 

and <l>i = V+Xt, $3 = V_Xt, eqs. (IT6T]) and (IT62|) take the final form: 

2FX3P^X, = 2F(P+X3P+X, - i?,G3 + i+X^Gt) + X3(2)+X2 - G?), ([161]) 
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2FX^VlXt = 2F{V^X;V^Xt - R^Gt + ^-^3^3) + X^iV^X'f - Gl). ^m 

On the other hand, in the corollary [7] of section [3] we obtained the system of 3rd-order 
equations, which should be an integrated version of equations of theorem [10] and other higher- 
order equations. Here we show that indeed all equations of theorem [10] follow from the ones 
of corollary [7] 

Proof. One has the following equations for the derivatives of Gt,G^, Aj^ and A^: 

XsV+Gt = V+X;Gt - RtV.Xt + UXzV+Xt, mM 

XsV.Gs = V.X.Gs - RiV+Xt + ^.X^V.Xt. ([225]) 

2FV_Gt = -(G3 - 2S,.F)V+Xt + {Gt + 2^+F)P_X„ (230) 

2FP+G3 = (G3 + 2i_F)V+Xt - [Gt - 2^+F)V_Xt, (231) 

XsV+A+ = V+X;A+ + 2a^{RtV_A^ - ^^X^V+A^), (232) 

XaV^A^ = V^X-iA^ + 2{R^V+A^ - i^Xj,V^A^). (233) 

2A^V.A+ = -(A_ + A^^A'^)V+A'^ + (A+ - A^+a^A'^)V^A'^, (234) 

2A^V+A_ = {A_ - Ai_A^)V+A^ - (A+ + Ai+a^A^)V_A^. (235) 

(Writing them out, we used e.g. formulas AFGa = GsV+Xt-GtV.Xt and AA^Ga = A_V+A^- 
Aj^V-A? .) Recall also relations: 

VlX^ = V+V^Xt + 2a'^V+V^A'^, V'^X^ = V+V^Xt + 2V+V^A^, 

V+V_X^ = V'iXt + 2a^V'iA^ = V^Xt + 2VlA^ , 
which are consequences of 

P+X3 = V^Xt + 2g'^V^A^, V^Xs = V+Xt + 2V+A'^, 

and which show that the written above are in fact the only five independent 4th-order equa- 
tions. The rest of the proof is long, tedious calculations, in which we take derivatives of 
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the 3rd-order PDE and compare them with such combinations of the 4th-order PDE that 
the highest derivative terms cancel out. Then we see that the rest reduces to the 3rd-order 
equations themselves. In this way one verifies e.g. that 



P+dni = V_Xt ■ ^62) + 2FX^V+Xt ■ ([ISO]), 
and, similarly, or by symmetry, 

P-dlZi = V+Xt ■ ^m) + 2FX^V_Xt ■ ^M). 
Differentiating (I174p . one can get a combination 

FX^V+^^ + cr'^V^A^ ■ ([ISID = 0, 
which turns out to be equivalent to 

a^ ■ ([T59D + X^V+A^ ■ ([163]) = 0. 
Then by symmetry one has also 

a^ ■ ^^ + X^V_A^ ■ ([US]) = 
from the combination 

FX^V^^^ + o'^V+A^ ■ ([162]) = 0. 

If one considers equation FX^V^iVt + V3) = and applies (l45ll and (I179p . one comes to the 
combination 

V+Xt ■ ([m]) + 2FX^V.Xt ■ ([ISO]) = 0, 
and analogously or by symmetry one also gets the combination 

V^Xt ■ ([ISD + 2FXsV+Xt ■ ([ISO]) = 0. 

Thus, one sees that all five independent 4th-order equations indeed follow from the 3rd- 
order equations. D 

References 

[1] M. Adler, T. Shiota, P. van Moerbeke. Random matrices, vertex operators and the 
Virasoro algebra. Physics Letters A, 208:67-78, 1995. 



50 



[2] M. Adler, P. van Moerbeke. Herinitian, symmetric and symplectic random ensembles: 
PDEs for the distribution of the spectrum. Ann. Math., 153:149-189, 2001. 

[3] M. Adler, P. van Moerbeke. The spectrum of coupled random matrices. Ann. Math., 
149:921-976, 1999. 

[4] M. Adler, P. van Moerbeke. PDEs for the joint distributions of the Dyson, Airy and Sine 
processes. Ann. Prob., 33:1326, 2005. 

[5] M. Bertola, B. Eynard, J. Harnad. Duality, biorthogonal polynomials and multimatrix 
models. Commun. Math. Phys., 229:73, 2002. 

[6] F. Dyson. A Brownian-Motion Model for the Eigenvalues of a Random Matrix. Journ. 
Math. Phys., 3:1191-1198, 1962. 

[7] B. Eynard, M.-L. Mehta. Matrices coupled in a chain: I. Eigenvalue correlations. Journ. 
Phys. A, 31:4449-4456, 1998. 

[8] P. Forrester. Log-gases and Random Matrices. www.ms.unimelb.edu.au/ mat- 
pjf/matpjf.html, online book, 2005. 

[9] Harish- Chandra. Differential operators on a semisimple Lie algebra. Amer. Journal 
Math., 79:87-120, 1957. 

[10] C Itzykson, J.-B. Zuber. The planar approximation II. Journ. Math. Phys., 21:411-421, 
1980. 

[11] M. Kardar, G. Parisi, Y.-C Zhang. Dynamic scaling of growing interfaces. Phys. Rev. 
Lett, 56:889-892, 1986. 

[12] T. Liggett. Interacting Particle Systems. [Reprint of the 1985 original.] Berlin: Springer- 
Verlag, 2005. 

[13] A. F. Macedo- Junior, A. M. S. Macedo. Brownian-motion ensembles: correlation func- 
tions of determinantal processes. Journ. Physics A, 41:015004, 2008. 

[14] M.-L. Mehta. Random matrices, 3rd edition. Elsevier, San Diego, 2004. 

[15] A. Newell. Solitons in Mathematics and Physics. Cambridge University Press, 1985. 

[16] M. Prahofer, H. Spohn. Scale invariance of the PNG droplet and the Airy process. Journ. 
Stat. Phys., 108:1071, 2002. 



51 



[17] I. Rumanov. The correspondence between Tracy- Widom (TW) and Adler-Shiota-van 
Moerbeke (ASvM) approaches in random matrix theory: the Gaussian case. Journal of 
Mathematical Physics, 49:043503, 2008. 

[18] I. Rumanov. Universal structure and universal PDE for unitary ensembles. Journal of 
Mathematical Physics, 51:083512, 2010. 

[19] I. Rumanov. All the lowest order PDE for spectral gaps of Gaussian matrices. 
arXiv:/1008.3560, 2010. 

[20] C.A. Tracy, H. Widom. Level-Spacing Distributions and the Airy Kernel. Commun. 
Math. Phys., 159:151-174, 1994. 

[21] C.A. Tracy, H. Widom. Fredholm determinants, differential equations and matrix models. 
Commun. Math. Phys., 163:33-72, 1994. 

[22] C.A. Tracy, H. Widom. A system of differential equations for the Airy process. Electronic 
Comm. Probab., 8:93-98, 2003. 

[23] C.A. Tracy, H. Widom. Differential equations for Dyson processes. Commun. Math. 
Phys., 252:7-41, 2004. 

[24] C.A. Tracy, H. Widom. Asymptotics in ASEP with step initial condition. Commun. 
Math. Phys., 290:129-154, 2009. 

[25] K. Ueno, K. Takasaki. The Toda Lattice Hierarchy. Adv. Stud. Pure Math., 4:1-95, 1984. 



52 



